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HIGHLIGHTS 


•  Simultaneous  input  and  state  estimation  proposed  for  solid  oxide  fuel  cell. 

•  The  stability  of  the  adaptive  observer  is  proven  using  the  Lyapunov  function  method. 

•  The  proposed  observer  is  robust  to  fluctuations  in  the  input  variables. 
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An  adaptive  nonlinear  observer  design  for  the  planar  solid  oxide  fuel  cell  (SOFC)  is  presented  in  this 
work.  This  observer  is  based  on  a  lumped  parameter  model  of  the  SOFC  and  it  can  simultaneously  es¬ 
timate  the  inputs  and  the  states  of  the  system.  Considering  the  inputs  as  unknown  parameters  is  ad¬ 
vantageous  because  some  of  the  input  parameters  are  not  practically  measurable  in  a  SOFC  stack.  The 
asymptotic  stability  of  the  proposed  observer  is  proven  using  the  Lyapunov  function  method  and  is  based 
on  the  concept  of  input-to-state  stability  for  cascaded  systems.  The  simulations  show  that  the  developed 
observer  can  track  the  temperature  and  species  concentration  profiles  in  the  planar  SOFC  during  step 
changes  in  the  input  variables  and  can  simultaneously  predict  the  input  variables.  The  adaptive  observer 
presented  is  valid  for  a  wide  operating  range,  requires  fewer  variables  to  be  measured,  and  is  robust  to 
fluctuations  in  the  input  variables. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  challenges  in  commercializing  solid  oxide  fuel  cells  (SOFC’s) 
may  be  classified  as  those  pertaining  to  material  science,  reaction 
kinetics  and  control  and  operation.  Even  if  the  materials  and  ki¬ 
netics  related  issues  are  sorted  out,  there  are  several  technological 
challenges  relating  to  control  and  operation  that  must  be  addressed 
before  successful  commercialisation.  One  of  the  important  issues  in 
the  planar  configuration  SOFC  is  the  non-uniform  spatial  temper¬ 
ature  distribution  in  the  cell,  which  places  serious  limitations  on  its 
performance.  The  non-uniformity  of  the  temperature  in  the  cell 
will  lead  to  hot  spots  and  thermal  stresses  that  in  turn  increases  the 
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probability  of  failure  and  degradation  of  the  cell  [1,2],  These  issues 
are  especially  predominant  during  the  transient  operation  of  the 
SOFC.  The  thermal  gradients  are  especially  high  in  the  case  of 
methane  fuelled  SOFC  with  internal  reforming.  Thermal  gradients 
would  also  be  significant  in  SOFC  systems  intended  for  load 
following  and  frequent  on/off  applications.  This  thermal  manage¬ 
ment  problem  has  an  important  bearing  on  the  efficiency,  life  and 
reliability  of  the  cell.  Therefore,  thermal  management  is  essential 
not  only  to  prevent  the  damage  to  the  fuel  cell  and  thereby  maxi¬ 
mise  the  cell  life,  but  also  to  improve  its  efficiency  and  performance 
[3,4],  For  effective  thermal  management,  information  about  the 
temperature  distributions  inside  the  cell  is  required.  Since  this  is 
generally  not  measurable,  dynamic  estimation  is  one  option  for 
obtaining  this  information. 

The  need  to  minimise  the  temperature  gradients  in  the  cell  has 
been  realised  in  the  past  and  methods  have  been  presented  [2]  for 
designing  control  strategies  to  minimise  the  spatial  temperature 
variation  in  a  molten  carbonate  fuel  cell.  A  numerical  bifurcation 
analysis  of  the  SOFC  was  performed  in  Ref.  [5],  which  gave  guide¬ 
lines  to  avoid  high  over-temperatures  in  the  SOFC.  Refs.  [6,7] 
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proposed  state  estimators  based  on  Kalman  filter  that  can  estimate 
the  temperature  and  the  concentration  profiles  inside  the  molten 
carbonate  fuel  cell.  In  Ref.  [8],  an  observer  was  presented  for  esti¬ 
mating  the  spatial  temperature  profiles  of  a  single  cell  SOFC,  which 
is  based  on  a  linearised  model.  Observers  designed  using  lineari¬ 
zation  techniques  are  not  very  suitable  for  the  SOFC,  which  is  a 
highly  non-linear  system. 

The  limited  amount  of  research  in  the  open  literature  use  line¬ 
arisation  approaches  to  observer  design  and  consider  a  single  cell 
design.  Many  models  in  the  literature  consider  the  cell  current  as 
the  input  [6,7]  and  treat  other  inputs  to  the  system  like  the  cathodic 
and  anodic  flow  rates  and  temperatures  as  constant  parameters. 
Some  designs  consider  the  inlet  mole  flow  rates  as  known  inputs 
that  are  measured  [8],  which  is  quite  possible  if  a  single  cell  is 
considered  or  if  the  whole  stack  is  lumped  together  as  a  single  unit 
for  the  design.  However,  considering  the  whole  stack  as  single  unit 
may  give  an  unrealistic  picture  of  the  temperature  distribution 
within  the  stack  because  in  a  planar  SOFC  stack,  different  cells  will 
have  different  temperature  profiles  [9]  both  during  steady  state  and 
transient  operation.  If  on  the  other  hand,  if  it  is  desired  to  include 
many  cells  into  the  model  based  observer  design,  it  may  not 
possible  to  measure  the  inlet  flow  rates  or  the  current  for  each  cell. 
In  this  situation,  it  is  beneficial  to  consider  the  inputs  as  unknowns 
and  estimate  their  values  in  addition  to  the  state  variables. 

The  SOFC  system  being  an  inherently  nonlinear  one  requires  a 
nonlinear  observer  for  the  purpose  of  a  wider  range  of  applicability 
and  accuracy.  Unlike  for  linear  systems,  there  is  no  systematic 
procedure  for  designing  a  state  observer  for  a  given  nonlinear 
model.  Empirical  observers  and  converging  observers  are  the  two 
kinds  of  observers  available  for  non-linear  systems  [10],  Kalman 
filter  is  a  common  converging  observer  which  is  very  useful  in  the 
case  of  linear  systems.  However,  its  extension  to  the  nonlinear 
system,  namely  the  extended  Kalman  filter  is  an  empirical  observer 
based  on  some  approximation  of  the  nonlinear  system  or  approx¬ 
imation  of  a  theoretical  best  estimation.  It  is  based  on  a  lineariza¬ 
tion  of  the  nonlinear  system  in  a  neighbourhood  of  its  estimation. 
Hence,  the  extended  Kalman  filter  is  not  a  globally  converging 
observer.  On  the  other  hand  theoretically  converging  observers  are 
only  applicable  for  a  certain  class  of  non-linear  systems  and  are 
usually  very  difficult  to  design.  Some  nonlinear  converging 
observer  designs  such  as  sliding  mode  observers  suffer  from 
problems  like  chattering.  Therefore,  a  Lyapunov  function  based 
adaptive  observer  design  has  been  undertaken,  which  is  a 
converging  observer.  Lyapunov  function  based  designs  suffer  from 
the  problem  that  there  are  no  systematic  methods  for  finding  a 
suitable  Lyapunov  function  for  a  given  system.  However,  in  this 
work,  the  prudent  choice  of  adaptive  parameters  and  the  adapta¬ 
tion  laws  allows  the  definition  of  a  Lyapunov  function,  using  which 
the  global  stability  of  the  error  dynamics  is  proven. 

In  a  previous  work  by  the  authors  [11],  an  adaptive  nonlinear 
observer  for  the  planar  SOFC  was  presented.  This  paper  presents 
an  extension  of  the  previous  work,  in  which  an  unknown  input 
adaptive  non-linear  observer  for  the  single  cell  SOFC  is  presented. 
It  is  assumed  that  the  inputs  (the  cell  current  density,  the  anode 
and  the  cathode  inlet  flow  rates  and  the  inlet  temperature)  are  not 
available  for  measurement  and  are  to  be  estimated  in  addition  to 
the  system  states.  It  is  also  assumed  that  the  stack  inlet  mole 
fractions  (which  will  be  the  same  for  all  cells  in  the  stack)  and 
temperature  as  the  measured  variables  used  in  the  estimation 
scheme.  Similar  to  the  previous  work,  the  asymptotic  stability  of 
the  proposed  observer  is  proven  using  the  Lyapunov  function 
method  considering  the  system  as  a  cascaded  one.  The  SOFC 
model  used  for  the  observer  design  is  discussed  in  Section  2.  In 
Section  3,  the  observer  design  for  the  system  is  presented.  In 
Section  4  the  stability  analysis  is  presented,  following  which,  the 


simulation  results  are  presented  in  Section  5.  Section  6  concludes 
the  paper. 

2.  Modelling 

The  SOFC  basically  consists  of  three  thin  layers  that  constitute 
the  PEN  (Positive  electrode/Electrolyte/Negative  electrode)  namely 
the  cathode,  the  electrolyte  and  the  anode.  The  PEN  is  sandwiched 
between  metallic  plates  called  interconnects  having  channels  for 
gas  flow  which  also  serves  as  current  collectors.  Hydrogen  is  sup¬ 
plied  on  the  anode  side  and  oxygen  on  the  cathode  side  and  the 
electrochemical  reaction  within  the  cell  produces  water  vapour  on 
the  anode  side.  The  reaction  is  exothermic  and  results  in  the  evo¬ 
lution  of  heat.  Fig.  1  shows  the  basic  construction  and  operation  of 
the  SOFC.  The  readers  may  refer  to  Refs.  [12,13]  for  further  details 
on  the  construction  and  functioning  of  the  SOFC. 

A  lumped  parameter  model  is  used  to  capture  the  thermal  dis¬ 
tribution  in  the  cell.  The  cell  model  is  constructed  using  6  lumps  for 
which  the  mass  and  the  thermal  balances  are  written.  Increasing 
the  number  of  lumps  will  result  in  more  accurate  prediction  of  the 
temperature  distribution.  This  model  for  the  hydrogen  fed  SOFC 
was  presented  and  validated  in  Ref.  [11]  and  more  details  can  be 
found  there.  For  the  benefit  of  the  reader,  the  model  development  is 
outlined  in  the  Appendix. 

An  improvement  in  the  model  used  in  this  work  compared  to 
Ref.  [11  ]  is  that  the  heat  capacities  of  the  anode  and  cathode  gases 
are  considered  as  functions  of  species  mole  fractions  at  the 
particular  lump  as  given  in  Eqs.  (A17— A20).  This  makes  the 
observer  design  more  complicated  but  will  improve  the  accuracy  of 
the  model. 

In  the  next  section,  an  unknown  input  adaptive  observer  design 
is  proposed  for  this  system.  The  cell  current,  the  anode  and  cathode 
inlet  flows  and  the  temperature  the  inlet  streams  will  be  consid¬ 
ered  as  the  inputs  that  need  to  be  estimated  in  addition  to  the 
states.  These  variables  will  be  defined  as  the  adaptive  parameters 
for  which  adaptation  laws  will  be  formulated  in  such  a  way  so  as  to 
ensure  the  asymptotic  stability  of  the  error  dynamics. 

3.  Adaptive  observer  design 

In  this  section,  an  unknown  input  adaptive  observer  based  on 
Lyapunov  function  approach  is  presented  for  estimating  the  states 
of  the  SOFC.  It  is  assumed  that  the  concentrations  of  the  three 
species  at  the  inlet  of  the  cell  are  measured  and  they  are  same  as 
those  in  the  first  lump.  Similarly,  it  is  assumed  that  the  temper¬ 
ature  of  the  first  lump  is  measured.  The  proposed  observer 
equations  are, 


Ref.  [14]). 


1220 


P  Vijay  et  al  /  Journal  of  Power  Sources  248  (2014)  1218-1233 


Pmol^gas^p  =  Vu/a  ~  4  4  ~  V H2  ^  -  E  4  (J)  +  E  4o  (?)  ^  “  4  (J) 

Pmoi^gas^  =  4*4  +  41  ^  -  E  4  0 + E  4o  0  j  -  44  -  42  ^  -  E  4  0 + E  4o  0  j  -  4  0 


(1) 

(2) 


where,  i  =  1  for  Eq.  (1)  and  i  =  2,... ,6  for  Eq.  (2). 


Pmol^gas  ^2°  -  “  4o7A  -  y^Q  |  “  E  4  (7)  +  E  4o  ( } )  )  +  rH20  (j) 


j= 1  j=l 


Pmol^gas  ^2°  =  4(4  +  4,0  (  “  E  4  (J )  +  E  4o  0  J  “  ^HjO^A  ~  4o  (  “  E  4  0  +  E  4o  0  )  +  rH20  0  (4) 


where,  i  =  1  for  Eq.  (3)  and  i  =  2,... ,6  for  Eq.  (4). 


pamoi^asd^2=yo4c -44" 


(5) 


f’mol^gas  —  yo2  +yo2  ) 


-E4 


0  -4^-4  fc-E4 


(6) 


where,  i  =  1  for  Eq.  (5)  and  i  =  2,... ,6  for  Eq.  (6). 


psVscSp~sf  - (44  +j/h20cp2°) t  ^+(44+^44  4 _  (44+4o40)rsfA 


-  (44  +4o4°)0 (-E4 0 +E4o0  j  -  (44  +44)04 -  (44 +44)0  E4 (7)j 

■W) 


-M4-ciors+- 


(0.5l\  IR, 

i  k 


(7) 
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=  (>W  +3?koCp2°)fs~1?A  +  (y'H2  Cp2  +ykoC»20)fs_1  (-E4®  +E4o(T)) 

+  (jV£2  +y$2c"2)?<r'Fc  +  (y'o2<£2  +y%^2)t1^  ^-£40^  -  (54cp2  +yb2o<£2°)rs£A 

-  (jV°2  +yN2Cp2)rsfc  -  (jW2  +  yko&°) %  (  -  ±  ^  (?)  +  £  40)) 


-  (y'o2cp2  +yN2Cp2)fs  (  -  £40 )  +^-?(ci  -  + 


where,  i  =  1  for  Eq.  (7)  and  i  =  2,..., 6  for  Eq.  (8).  The  parameter  and  state  estimation  errors  are  classified  into  the  following 
adaptation  laws  are  given  as:  groups: 


T?  "  +  ~*'hE,0 

d?  ^3  ?  c  (^O  -  C  yH2S 

- ^ 

fjt  c  (vh2o  -ph2)S  Ph2oS\ 

-  yH20£yi2Q — 2f — + %<rxr) 

i  1  /V  F  v02l_F  v(hs\ 

«as  Vy°2  A  2F  4  2fJ 

(«  +y&0tf% 

+  Fc  (yo2Cp2  +yN2cp2)£rs1) 


(9) 

(10) 


ep-  %  %0  %  £/]T,  (13) 

E2  =  | 'EyH2  EyH20  £y^]T,  where  i  =  2,.. .,6.  (14) 

E3  -  [£T>  Ejin  j  ^  (15) 

E4  m  [ev?T.  where  i  m  2,  ...,6.  (16) 

where  Ep„ ,  Ep„ ,  Ey,  ,  £y,  ,  Eyi  ,  Ep ,  £/,  Er,„  are  the  parameter  and 
the  state "estimatiori  errors  giv^n  asS; 

£Fjn  =  F’"  -  Fa  (17) 

£f?,=F“-Fc  (18) 

£yH2  =  y[r2  -  ?h2  (19) 


where  Ai,  A2,  A3,  A4  >  O.TheA's  in  the  adaptation  laws  are  the  tuning 
parameters  introduced  to  influence  the  convergence  time  for 
the  estimation,  since  only  the  asymptotic  stability  of  the  error 
dynamics  is  assured.  Therefore,  the  tuning  parameters  will  help 
to  achieve  faster  settling  time.  This  can  be  understood  as 
being  analogous  to  decreasing  the  capacitance  of  a  thermal  or 
electrical  system,  thereby  achieving  a  decreased  time  constant. 
The  adaptation  laws  are  designed  in  such  a  way  that  the  error  dy¬ 
namics  become  asymptotically  stable  as  described  in  the  next 
section. 

4.  Stability  of  the  error  dynamics 

In  this  section,  the  stability  analysis  for  the  error  dynamics  is 
presented  using  the  Lyapunov  function  method.  For  the  purpose 
of  proving  the  asymptotic  stability  of  the  system,  the  parameter 


£yH2o  =  y^o  -  y'tho 

(20) 

%  =  yo2  -  yb2 

(21) 

En  =  T’  -  ?s 

(22) 

£,  =  /-/ 

(23) 

• 

Epn  =  Tm-T 

(24) 

From  model  equations  in  Eqs.  (A21)— (A28)  and  observer  equa¬ 
tions  in  Eqs.  (1)— (12),  the  error  dynamics  are  derived  as  follows: 


(25) 


iX,"^  -  \ (ly  +  ^-^d- iy)  +y<„,’  **  +3fe'£^5k)?(,-  „ . (fjf  +  (>H,o-^".)S(0f) 

where,  i  =  1  for  Eq.  (25)  and  i  =  2 . 6  for  Eq.  (26). 

Pmol^gas  df2°  =  £f?yi?20  -  £yi,20  (fA  +  ^  °2f  ~~  ^l)  ~  ^H»2oEF-  ~  5^0^  **°2F  ^  ^  +  ^  2f^  ®  (27) 

pU^as^  =  £y'2'0(Fr+Ko2Fi'H2)s  a-D/)  +&* + (f  - 1)  -  ^  (Er + 

-  ykoEtf  -  ykoE<('VH2°2Fv"2)S  ® + E'{VH2F]S  (i)  (28) 

where,  i  =  1  for  Eq.  (27)  and  i  =  2 . 6  for  Eq.  (28). 
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pamorJE^  =  E^-^a- i)/)  +  y'2\„ -y^Ei^a-i)-^^-^©/)  -y^  (30) 

where,  i  =  1  for  Eq.  (29)  and  i  —  2 . 6  for  Eq.  (30). 
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(38) 

(35) 
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The  stability  analysis  for  the  proposed  observer  is  based  on  the 
concept  of  input-to-state  stability  for  cascaded  systems  [15],  The 
asymptotic  stability  of  the  system  of  equations  £5  is  proven  first. 
The  asymptotic  stability  of  unforced  system  E2  with  the  system  Ei 
as  the  input  is  then  proven.  Then,  considering  systems  Ei  and  E2  as 


E3  =/2(E1,E2,E3) 

(39) 

e4  =/2(e,,e2,e3,e4) 

(40) 

The  global  uniform  asymptotic  stability  of  the  origin  of  the 
cascaded  system  described  in  Eqs.  (37)  and  (38)  is  assured  [15]  if, 
the  origin  of  Eq.  (37)  is  globally  uniformly  asymptotically  stable  and 
the  origin  of  the  Eq.  (38)  with  Ei  as  input  is  input-to-state  stable.  If 
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system  comprising  of  Eq.  (39)  is  proved  to  be  input-to-state  stable 
with  Ei  and  E2  (which  are  themselves  asymptotically  stable)  as 
inputs,  then  the  system  comprising  of  Eqs.  (37)-(39)  is  globally 
uniformly  asymptotically  stable.  Similarly,  if  the  asymptotically 
stable  system  comprising  of  Eqs.  (37)— (39)  are  considered  as  inputs 
to  the  system  in  Eq.  (40)  and  the  origin  of  Eq.  (40)  is  proved  to  be 
input-to-state  stable,  this  means  that  the  complete  system  of 
equation  from  Eq.  (37)-(40)  is  globally  uniformly  asymptotically 
stable. 


adaptation  laws  for  the  three  adapted  variables  (viz.  FjJ1,  FJ P,  I)  is 
such  that  the  terms  involving  the  corresponding  error  terms  in  V 
are  cancelled.  It  is  apparent  that  the  derivative  of  the  Lyapunov 
function  is  negative  definite  since  the  terms  inside  the  brackets  are 
positive  quantities.  This  proves  the  global  uniform  asymptotic 
stability  of  the  origin  of  the  set  of  equations  in  Eq.  (37). 

Next,  Eq.  (38)  is  examined  as  an  unforced  system  where  the 
inputs  from  Eq.  (37)  are  zero  in  order  to  establish  its  input-to-state 
stability.  Eliminating  the  terms  corresponding  to  Ei  in  Eq.  (38),  a  set 
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To  prove  the  asymptotic  stability  of  the  system  comprising 
Eq.  (37),  consider  the  following  Lyapunov  function  candidate: 


It  can  be  seen  that  V  =  0,  when  the  errors  are  zero  and  positive 
otherwise.  Therefore  the  Lyapunov  function  is  positive  definite.  The 
function  is  also  continuously  differentiable  and  radially 
unbounded. 

The  derivative  of  the  Lyapunov  function  is  given  as: 


'  ' —  >yr-  7  F pin  “f“  „  Erin  -f-  Ej  “f~  Eyl 

3£Fi„  OFfin  rc  0£ ,  dEy ,  y»2 


0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

-a5  0  0  0 

a5  -a6  0  0 

0  0  -b2  0 

0  0  b2  -b3 

0  0  0  b3 

0  0  0  0 

0  0  0  0 


(44) 


of  linear  differential  equations  is  obtained  that  can  be  written  in 
matrix  vector  notation  as  follows: 

where, 


“2  =  Fa  -E  4 (0  +  E  4oC0  (45) 

j=l  J=1 

°3  =  4-E4(o+E4q(o  (46) 

a4  =  fa  -  E  4 (o + E  4o(0  (47) 


v  =  ~2Ek  (fa -E 4(0 + E 4o(0  -  2Ek  0  \fa  -  E 4(0 + E 4o(o  -  (  f-  -  e 4(o 


(43) 


which  after  substitution  becomes: 

It  can  be  noted  that  the  derivative  of  the  Lyapunov  function  in 
the  form  in  Eq.  (43)  is  obtained  because  the  formulation  of  the 


“5  =  fr-E4(0+E4o(0 

j=t  j=i 


(48) 
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a6  =  Fin-E'JH2(/)  +  E4o« 


,dEr  .in 

Psvs4sr  =  +fa£7 

■i"(yH2Cp2  +yi?20<£2°)  +  (yoVp2  +yw2cp2)^ 

+  Th2o 

2 

h = E-E^w 

j=\ 

(50) 

=  ^n-E4(/) 

j=  1 

(51) 

b4  =  E-E^w 

j= 1 

(52) 

=  E-EE2(/) 

(53) 

=  E-ErJo2(/) 

(54) 

The  constants  ai . be  are  composed  of  known  quantities  and 

it  can  be  verified  that  the  matrix  in  Eq.  (44)  has  real  negative 
eigen-values.  Therefore,  the  unforced  system  described  in  Eq. 
(38)  is  globally  exponentially  stable.  This  proves  the  input-to- 
state  stability  of  the  origin  of  the  system  in  Eq.  (38)  and  the 
global  asymptotic  stability  of  the  system  comprising  of  Eqs.  (37) 
and  (38). 

In  the  following,  the  input-to-state  stability  of  the  system  in  Eq. 

(39)  is  proven  considering  Ei  and  E2  as  inputs  to  it.  Eliminating  the 
terms  corresponding  to  Ei  and  E2  in  Eq.  (39),  the  following  differ¬ 
ential  equation  is  obtained  for  E3. 

The  following  Lyapunov  function  is  considered  for  the  stability 


-  Er;  (yW-F?  +yW°f?  +yM-(l'",V,,)s' 

(55) 


(56) 


It  can  be  seen  that  V  =  0,  when  the  errors  are  zero  and  positive 
otherwise.  Therefore  the  Lyapunov  function  is  positive  definite.  The 
function  is  also  continuously  differentiable  and  radially  unbounded. 

The  derivative  of  the  Lyapunov  function  is  given  as: 

l/  =  -2 

+  ^o^0(^  +  K°2f^)S/)  (57) 

+J'o2<£2  (f'c  +  +yi?2Cp2  (Fc  +  5§r  ')) 

It  can  be  noted  that  the  derivative  of  the  Lyapunov  function  in  the 
form  in  Eq.  ( 57 )  is  obtained  because  the  formulation  of  the  adaptation 
laws  for  the  adapted  variables  tn  is  such  that  the  terms  involving  the 
corresponding  error  terms  in  V  are  cancelled.  It  is  apparent  that  the 
derivative  of  the  Lyapunov  function  is  negative  definite  since  the  term 
inside  the  bracket  is  clearly  a  positive  quantity.  This  proves  the  input- 
to-state  stability  of  the  origin  of  the  system  in  Eq.  (37)  and  the  global 
asymptotic  stability  of  the  system  comprising  of  Eqs.  (37)-(39). 

Finally,  in  the  following,  it  is  proven  that  the  system  E4  is  input- 
to-state  stable  with  systems  Ej,  E2  and  E3  considered  as  inputs  to  it. 
Imposing  this  condition,  the  following  linear  differential  equations 
are  obtained, 
where, 


— c5  +  d5 


C2  - yiW (f?  +  '"l)s')  +yH,oCp’0(f»  *' +  &<? (f? - ^1)  +yfe4f  (f?  + 


(59) 
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All  the  diagonal  terms  in  the  matrix  are  negative  and  other 
terms  are  positive.  It  can  be  verified  that  the  matrix  in  Eq.  (58)  has 
real  negative  eigen-values,  making  the  unforced  system  described 
in  Eq.  (40)  is  globally  exponentially  stable.  This  proves  the  input-to- 
state  stability  of  the  origin  of  the  system  in  Eq.  (40)  and  the  global 
asymptotic  stability  of  the  system  comprising  of  Eqs.  (37)— (40). 
Note  that  the  observer  design  elaborated  in  this  work  is  applicable 
for  models  with  any  number  of  lumps. 


5.  Simulation  results  and  discussion 

In  this  section,  the  simulated  responses  for  the  model  and  the 
observer  states  subjected  to  step  changes  in  the  inputs,  viz.,  the 
current  density,  the  inlet  temperature  and  the  two  inlet  flow  rates, 
are  presented.  In  addition,  the  response  of  the  observer  inputs  that 
are  considered  as  the  adaptation  variables,  are  also  presented. 
These  responses  correspond  to  the  step  change  in  the  input  vari¬ 
ables  at  500  s  as  follows:  step  change  in  inlet  temperature  from 
1000  K  to  1100  K;  step  change  in  the  current  density  from 
600  A  m~2— 550  A  m-2;  step  change  in  anode  inlet  flow  from 
0.0001  mol  s-1  to  0.0002  mol  s-1 ;  step  change  in  cathode  inlet  flow 
from  2.742e~3  mol  s  1  to  3e  3  mol  s_1.  The  fuel  flow  rate  of 
0.0001  mol  s  1  corresponds  to  a  value  of  400  ml  min  1  and  is 
adopted  from  Ref.  [16],  A  high  step  increase  in  the  fuel  flow  rate  is 
considered  so  as  to  demonstrate  the  relevance  of  the  observer  over 
a  wide  operating  range  and  at  low  fuel  utilisations.  The  fuel  flow 
rate  is  especially  susceptible  to  much  fluctuation  if  a  reformer  is 
involved  in  the  upstream.  The  operating  conditions  and  geomet¬ 
rical  parameters  are  taken  from  Ref.  [16],  Some  of  the  important 
parameters  used  in  the  simulations  are  given  in  Table  1. 

The  responses  of  the  model  states  are  represented  in  Figs.  2-5  in 
full  lines  and  the  responses  of  the  observer  states  are  represented  in 
dashed  lines.  Only  the  results  from  nodes  2, 4  and  6  are  presented 
for  clarity.  It  can  be  seen  that  the  observer  is  able  to  track  the  state 
variables  across  all  the  six  nodes.  In  Fig.  3,  the  species  mole  frac¬ 
tions  are  not  changing  after  the  step  change  because  the  change  in 
inlet  temperature  has  no  effect  on  the  species  mole  fractions. 
Similarly,  as  seen  in  Fig.  4,  the  change  in  anode  inlet  flow  has  no 
effect  on  the  oxygen  mole  fraction.  From  Fig.  5,  it  is  also  evident  that 
the  change  in  cathode  inlet  flow  has  no  effect  on  the  species  mole 
fractions  on  the  anode  side. 


Values  of  some  important  parameters  used  in  the  simulation. 


Parameter 

Value 

Anode  catalyst  width,  m 

0.25  x  to  3 

Cathode  catalyst  width,  m 

0.03  x  10  3 

Porosity 

0.4 

Activation  energy,  J  mol  1 

120,000 

Channel  height,  m 

0.75  x  10“3 

5  x  10“3 

Electrolyte  width,  m 

1  x  10“5 

Ohmic  resistance,  £2  m2 

2.94  x  10-6 

Limiting  current  density,  A  m  2 

3500 

PsCp  Q  m-3  K-1) 

106 
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Fig.  2.  Observer  tracking  during  step  change  in  cu 


As  discussed  earlier,  some  tuning  variables  (A’s)  were  introduced 
so  as  to  achieve  faster  settling  of  the  error  variables.  In  general,  the 
higher  the  value  of  the  tuning  variable,  the  faster  is  the  convergence 
of  the  adaptation  variables.  The  simulations  revealed  that  any  one 
the  tuning  parameter  also  had  an  influence  on  the  convergence  of 
other  variables  in  the  system.  In  order  to  understand  this  let  us 
consider  the  following  error  dynamics  equation: 
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Fig.  4.  Observer  tracking  during  step  change  in  anode  inlet  flow. 


Fig.  5.  Observer  tracking  during  step  change  in  cathode  inlet  flow. 
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Following  the  method  in  Ref.  [17],  Eq.  (69)  c; 
follows. 
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where  P  is  the  Laplace  operator.  Let  us  assume  that  Ey i  and  £;  are 
inputs  to  Eq.  (70).  Taking  into  consideration  that  Eq.  (7&)  should  be 
equal  to  zero  at  steady  state  and  assuming  that  only  Ej  has  not 
settled  to  zero,  the  equation  reduces  to  the  following. 


-  >1^ 


(71) 
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It  is  evident  from  Eq.  (71),  that  under  the  given  circumstances, 
it  is  possible  that  EF,„  has  a  non-zero  solution  and  it  can  settle  into 
this  solution  at  steady  state  instead  of  the  origin.  Similar  sort  of 
arguments  can  be  developed  for  other  variables.  However,  if  the 
inputs  are  in  the  form  of  a  persistent  excitation,  parameter 
convergence  can  be  guaranteed  [17],  But  in  the  fuel  cell,  the  inputs 
are  usually  constant  with  possible  step  or  ramp  changes.  Therefore, 
it  can  be  understood  that  the  coupled  nature  of  the  error  dynamics 
makes  it  sensitive  to  the  tuning  parameters.  In  view  of  this,  the 
tuning  parameters  were  chosen  as  large  as  possible.  Satisfactory 
convergence  was  obtained  using  the  values  for  the  tuning  variables 
as: A3  =  le25;  Ai  =  A2  =  A4  =  le4.  Any  larger  values  resulted  in 
problems  in  the  numerical  simulations.  These  large  values  for  the 
tuning  variables  increased  the  simulation  time.  The  simulation  of 
the  model  and  the  observer  for  1000  s  using  Matlab  Simulink 
(Version  7.8.0  R2009a)  in  a  personal  computer  with  Intel  Core  2 
Duo  processor  with  a  speed  of  3  GHz  and  1.93  GB  of  RAM  took 
about  836  s  in  real  time.  Therefore,  the  proposed  observer  can  be 
implemented  in  real  time. 

The  quantities  plotted  in  Figs.  6—9  are  the  estimated  values  of 
the  input  variables  and  are  not  the  actual  physical  values  of  these 
variables.  As  can  be  seen  form  Figs.  6-9,  some  adapted  input  var¬ 
iables  have  a  small  deviation  when  step  change  is  effected  to  other 
input.  For  example,  in  Fig.  6,  the  step  change  in  current  density 
causes  a  small  deviation  in  the  estimated  value  of  cathode  inlet 
flow  (the  error  is  in  the  order  of  2e-9).  Similarly  there  is  an  error  in 
the  anode  inlet  flow  of  the  order  of  9e-4.  Again  in  Fig.  7,  there  are 
small  errors  in  the  anode  and  cathode  inlet  flows  of  the  orders  of 
le-13  and  le-14,  respectively.  Also  note  that  the  errors  in  estimated 
inputs  could  be  either  positive  or  negative  (for  example,  in  Fig.  6, 
the  error  in  the  anode  and  cathode  inlet  flows  is  positive,  while 
they  are  negative  in  Fig.  7).  However,  these  deviations  are  too  small 
to  have  any  considerable  effect  on  the  state  variable  estimation  as 
seen  in  the  simulations. 

As  the  tuning  parameters  influence  the  convergence  of  these 
estimated  variables,  large  values  for  these  parameters  reduce  the 
error  in  the  input  variables.  However,  this  can  also  cause  spikes  in 
the  estimated  response  as  seen  in  Fig.  8.  This  seems  to  be  a 
disadvantage  of  considering  more  adopted  variables  because  less 
number  of  these  variables  did  not  result  in  any  error  [11],  Smaller 
values  for  the  tuning  parameters  reduce  the  overshoot,  but  increase 
the  error  in  the  input  estimation. 

6.  Conclusions 


parameters  in  the  observer  design.  This  method  also  offers  scope 
for  extending  the  observer  for  the  estimation  of  distributions  in  the 
entire  stack.  The  non-linear  model  based  observer  presented  in  this 
work  is  also  robust  to  the  disturbances  in  inputs  and  measurement 
noises.  These  input  disturbances  can  occur  in  the  system  because  of 
the  performance  issues  with  the  balance  of  plant  components  (like 
compressors,  power  electronics  and  heat  exchangers)  in  the  circuit. 
This  approach  will  be  extended  in  a  future  work  for  estimation  of 
temperatures  in  a  multi  cell  planar  SOFC  stack. 

Appendix 

The  model  development  for  the  purpose  of  observer  design  is 
outlined  in  this  Appendix.  The  hydrogen  balance  for  the  anode 
channel  of  the  SOFC  results  in  the  following  set  of  equations. 


pLK^  =  yW  -  yUA  -  4  (0  (ai  ) 

=  4Fa  1  -  VhA  -  4  (0  (A2) 

where,  i  =  1  for  Eq.  (Al)  and  i  =  2,. ..,6  for  Eq.  (A2). 

Similarly,  the  water  vapour  balance  for  the  anode  gas  channel  is 
given  by  the  following  set  of  equations. 

Pmol^gas^  =  TuVa  -  4 o4  +  4o(0  (A3) 

Pin ol4s^jr  =  4oFA  1  ~  4oFi  +  4o CO  (*4) 

where,  i  =  1  for  Eq.  (A3)  and  i  =  2,. ..,6  for  Eq.  (A4). 

The  oxygen  balance  in  the  cathode  channels  of  the  SOFC  are 
written  as  follows: 

p™i4sddt2  =  yoA  -  ybA  -  4  (o  (as) 

=  y'oAc  '  -  ybA  -  4  (0  (A6) 


In  this  paper,  an  adaptive  nonlinear  estimation  scheme  for  the 
simultaneous  estimation  of  states  and  inputs  in  a  planar  SOFC  is 
presented.  The  observer  design  is  based  on  an  improved  lumped 
parameter  model  of  the  SOFC  that  is  capable  of  providing  the 
spatial  profiles  of  temperatures  and  species  concentrations.  In  this 
design,  the  inputs  such  as  the  anode  and  cathode  flow  rates,  the 
inlet  stream  temperature  and  the  cell  current  are  considered  to  be 
unknown.  These  inputs  variables  are  estimated  in  addition  to  the 
state  variables.  The  stability  of  the  proposed  observer  is  proven 
using  the  Lyapunov  function  method  and  is  based  on  the  concept  of 
input-to-state  stability  for  cascaded  systems.  The  simulations  show 
that  the  developed  observer  can  track  the  temperature  and  species 
concentration  profiles  in  the  planar  SOFC  during  step  changes  in 
the  input  variables  and  can  simultaneously  estimate  the  inputs.  The 
observer  design  methodology  presented  in  this  paper  is  applicable 
for  models  with  any  number  of  lumps. 

The  unknown  input  non-linear  adaptive  observer  developed  in 
this  work  requires  fewer  variables  to  be  measured,  is  valid  for  a 
wide  range  of  operating  conditions  and  is  robust  to  variation  in  the 
input  variables  because  these  have  been  considered  as  adaptive 


where,  i  =  1  for  Eq.  (A5)  and  i  =  2,.  ..,6  for  Eq.  (A6). 

The  thermal  balances  for  the  solid  volume  of  the  SOFC  are  given 
in  the  following  equations: 


PsVs^  A  =  Fafa  +  W  -  4fa  -  A'c  +  AHr4(7)  -  V‘cenI 

(A7) 


h'A  'F{  1  +  1  -  hAFA  -  h‘cF[  +  AHr4(0 

~  Kell1 

(A8) 


where,  i  =  1  for  Eq.  (A7)  and  i  =  2,. ..,6  for  Eq.  (A8). 

The  reaction  rates  in  the  preceding  equations  are  given  as: 


4  CO  = 


Vh2I 


(A9) 
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r'ojn  =  2°/  (A10) 

4oW  =  Vj^r  (A11) 

The  flows  in  the  anode  and  cathode  channels  are  given  as 
follows: 


fi  =  fin-Ei(J)+E4oW 
4  =  Fc-^oJI) 


pL.E^=44n-44n-4  ( -£4«+£4o<0 


i=i  i=i 


-4(0 


pLK*sdft2  =  ^'fa  +  41  (  -  E 4(0 + E 4o(0 


-  th/a  -  4  y  -  E 4(0 + E  4o(0 

-4(0 


The  cell  voltages  in  each  of  the  lumps  are  given  by, 


Ell  =  Ecv  -  'Rohm  -  -  ^-h  1  (0/o5/) 


where  the  open  circuit  voltage  is  given  by  the  Nernst  equation  as: 


vocv 


4  + 


■  /  •  \0.5 

(ykM_\ 
V  4o  ) 


(A15) 


In  Eq.  (A15),  Eq  is  the  standard  cell  potential  given  as  [19]: 


4  =  Cl  -  ClT's  (A16) 

The  second,  third  and  the  fourth  terms  on  the  right  hand  side  of 
Eq.  (A14)  represent  the  ohmic,  concentration  and  the  activation 
over-potentials  [18]  that  are  the  sources  of  potential  losses  in  the 
fuel  cell.  The  simplified  expression  for  the  concentration  over¬ 
potential  involves  a  limiting  current  and  is  taken  from  Refs. 
[19,20],  The  expression  for  the  activation  over-potential  is  obtained 
from  the  Butler-Volmer  equation  if  the  transfer  co-efficients  (a) 
are  0.5,  which  is  usually  the  case  [19,20],  The  ohmic  resistance 
(Kohm)  and  the  exchange  current  (/o)  are  considered  constants  in 
this  work. 

The  specific  enthalpies  in  Eqs.  (A7)  and  (A8)  are  given  by, 


hA  =  4T’in  = 

(44  +4o4°)7'in 

(A17) 

4  =  <Es  = 

(Th242  +yH20Cp2°)TS 

(A18) 

hm  =  cCTin  = 

(44  +y[?2Cp2)7'in 

(A19) 

4  =  44  *■'! 

(44  +44)7s 

(A20) 

Substitution  of  Eqs.  (A9)-(A20)  in  Eqs.  (A1)-(A8)  results  in  the 
following  set  of  equations  for  the  SOFC. 


where,  i  =  1  for  Eq.  (A21 )  and  i  —  2,... ,6  for  Eq.  (A22). 

Similarly,  the  water  vapour  balance  for  the  anode  gas  channel  is 
given  by  the  following  set  of  equations. 


PmolEs 


d4o 

dt 


4o4"-4o41 

-  4o  (  -  E  4(*> + E  4oco)  +  4oC0 

(A23) 


Emol^gas 


d4( 

dt 


4o4"  +4o 


-£4(0+£4o(0 


-4o4n-4o 

+  rH2o(0 


-£4(0  +  £4o(0 

J=1  1=1  y 


(A24) 


where,  i  =  1  for  Eq.  (A23)  and  i  =  2,... ,6  for  Eq.  (A24). 

The  oxygen  balance  in  the  cathode  channels  of  the  SOFC  is  given 
as: 


pamoiEs-ar  =  EE  -  44n  - . 4  ^  -  E  4 (Oj  -  4 (0 

(A25) 

pamoiEsddt2  =  41E+41(-E4m)  -44" 

-4(  -E4(J)|  -4(0 


where,  i  =  1  for  Eq.  (A25)  and  i  =  2,... ,6  for  Eq.  (A26). 

The  thermal  balances  for  the  solid  volume  of  the  SOFC  are: 
where,  i  =  1  for  Eq.  (A27)  and  i  =  2,. ..,6  for  Eq.  (A28).  The  Eqs. 


+  (y$,$  +yi^-)T'“F?  -  (*W‘  +C-r)’» 

-  (j4,cS-+/„,0cr)i’s  ( -t  4w + E  Cm)  -  (44-  +<<$‘)rH? 

-  (y0,c?  +Acf)n(- 1  Co)  +^-i(4-£S,U+5i"(^-)) 

os^scff  -  (/rf  +C^°)T’f?+ (A,#  +y„,oc?0)is-' (-E4m +ECo>) 

+  (A,<?  +&#)■*'<?+ (44-  +<<?)*'  (-Eim)  -  (A.4-  ♦CxT)^ 

-  (y,,^  +<^)4f?  -  (yH.cs- +y„.„cr)7i  ( -  ±  Co + tcm) 

-  (A.CS-  +yW)n(-  ECo)  +  -  ci^^M-)) 


(A27) 


(A28) 


(A21)— (A28)  represent  the  SOFC 1 
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